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Abstract 

We show that there is a one-to-one correspondence between wave 
functions and surfaces in the position-momentum phase plane bounded 
by a closed curve satisfying an exact quantum condition refining the 
usual EBK condition. This is achieved using an old forgotten idea of 
Enrico Fermi. 

1 Introduction 

We address in this paper the following simple questions: 

Given a wave function ip(x) defined on the real line, is it 
possible to give an unambiguous two-dimensional pictorial rep- 
resentation of that function as a surface in position-momentum 
phase plane? Conversely, under which conditions can one asso- 
ciate to such a surface a wave function? 

We will show that there is indeed a one-to-one correspondence between 
wave functions and surfaces in phase plane whose boundary is a closed curve 
satisfying a certain quantum condition. Our proof is based on two results, 
the first of which has almost sunk into oblivion: "Fermi's trick" [5j which 
associates to every wave function a curve, and the existence of an exact 
quantization condition for such curves. More precisely, we will see that: 
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• Every twice differentiable complex function ip satisfies a linear (second 
order) differential equation 




f{x)\ + g{x) ^ = 



(1) 



• Closed phase plane curves 



P = P + {x) > , p = p (x) < {xa < x < xb) 



(2) 



are in one-to-one correspondence with functions ip provided they sat- 
isfy an exact quantization condition [T71 Q2J [2D] . 

Acknowledgements. The present work has been supported by the grant 
P20442-N13 of the Austrian Research Agency FWF. 

2 Fermi's trick 

In a largely forgotten paper [5] from 1930 Enrico Fermi shows that one 
could associate to every quantum state tp a certain curve gp(x,p) = in 
phase plane. He remarked that one could explicit determine a second or- 
der differential equation satisfied by a twice differentiable complex function. 
Fermi's work has recently rediscovered by Benenti and Strini [H [3]. The 
underlying idea is actually surprisingly simple. It consists in observing that 
any complex twice continuously differentiable function xfj(x) = R(x)e lS ( x ^ h 
(R(x) > and S(x) real) defined on configuration space trivially satisfies 
the partial differential equation 



(it is assumed throughout that R satisfies the concavity condition R" < 0). 
Performing the gauge transformation 




(3) 




d 



d 



this equation is equivalent to 
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(the equivalence of Eqns. ([3]) and (jH) can also be verified by a direct explicit 
calculation). The operator 

®-(-"s-*) >+ "if (5) 

appearing in the left-hand side of Eqn. (J3|) is the quantization of the real 
function 

9F(x,p) = { P -S') 2 + h 2 ^ (6) 

and the equation gp(x,p) = in general determines a curve Cf in the phase 
plane, which Fermi ultimately identifies with the state ip itself. Notice that 
the Fermi function only depends on the state tp in the sense that if the 
replacement of ip with Aip (A / a complex constant) changes neither gp 
nor g^. 

Let us illustrate Fermi's trick when ip is a Gaussian ip a ,b( x ) = e~( a+lb ^ x ' 2fi . 
In this case we have S(x) = —bx 2 /2 and R(x) = e~ ax2 / 2h hence 

gF(x,p) = (p + bx) 2 + ax 2 — ah. 

The Gaussian ip a ^ is thus a solution of the eigenvalue problem 

tp = —ahifj. 

For instance, if a = 1 and 6 = one recovers the fact that the standard 
Gaussian e~ x / 2h is the ground state of the harmonic oscillator. 



-ih— — h bx 
ax 



+ ax 



3 An exact quantization rule 

Consider the one-dimensional stationary Schrodinger equation 

h 2 d 2 



2m dx 2 



ifj(x) = [E - V(x)]i;(x). (7) 



We assume that the potential V is piecewise continuous, and that there exist 
exactly two real values xa and xb ("turning points") such that 

V(x) > E for — co < x < xa or xb < x < +oo 
V(x) = E for x = xa or x = xb 
V(x) < E for xa < x < xb- 
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It follows from elementary functional analysis (Sturm-Liouville theory) that 
the equation ([7]) has non-zero solutions only for a set of discrete values 
Eq < E\ < • • • of the energy E. In |17t [T8] Ma and Xu have shown that 
these values can be explicitly calculated using an exact quantization rule. 
The argument goes as follows: let x = V'VV' be the logarithmic derivative of 
the wavefunction ifj; a straightforward calculation shows that Schrodinger's 
equation ([7]) is equivalent to the Riccati equations 

-x'(x) = k(x) 2 + x(x) 2 for V(x) < E 
- X '(x) = -k(x) 2 +x(x) 2 for V(x)>E; 
where k{x) is given by 

-y / 2m(E - V(x)) for V(x) < E 

*(*) = < \ (8) 

-y/2m(V(x) - E) for V{x) > E 

(the function p(x) = hk(x) is the momentum ). Supposing that the potential 

V(x) is continuous at the turning points xa and xb, the exact quantum 

condition of Ma and Xu is then 

rx B j-XB 

k{x)dx - / x(x)[x'(x)}- 1 k'(x)dx = Ntt (9) 

where N is the number of nodes of x f° r V(x) < E (hence N—l is the number 
of nodes of ip in that region). Notice that if the second integral is neglected, 
the formula above reduces to the familiar approximate Bohr-Sommerfeld or 
WKB prescription 

rx B 

/ k{x)dx = Ntt. (10) 

Jx A 

We mention that Barclay [2] had shown in some older and unfortunately 
rather overlooked work that one can obtain exact quantization rules by a 
clever re-summation procedure of the higher order terms in the WKB series. 

A crucial fact noted by Qian and Dong [20] (also see Serrano et al. 
[2"H [22] ) is that condition Q can be put for all exactly solvable systems in 
the simple form 

XB fXB 

k(x)dx- k (x)dx = (N - l)vr (11) 

X A Jx A 

where the function /co is defined by ([8]) with E = Eq (the ground energy 
level) : 

hk (x) = y / 2m(E - V{x)) for V{x) < E. (12) 
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The exact quantization condition thus becomes 

r x B (-x B jyh 

/ y/2m[E - V{x)]dx = / y / 2m[E - V(x)]dx + — (13) 

Jx A Jx A * 

which shows that the energy levels E\ , E2, ■ ■ ■ are determined by the ground 
state energy level Eq. Denoting by C (resp. Co) the curve 

— p 2 + V(x) = E (resp. — p 2 + V(x) = Eq) 
2m 2m 

this is equivalent to 

/ pdx = I pdx + Nh. (14) 
JC Jco 

While a general condition for systems with n degrees of freedom is still 
lacking, the procedure outlined above still applies to n-dimensional systems 
with spherical symmetry in which case V(x) is replaced by the effective 
potential V e g(r) (see [TBJ [T71 [181 12Q] ) ■ This allows, in particular, to recover 
the energy levels of the hydrogen atom. 

4 Wave-functions and surfaces in phase plane 
4.1 From ^ to Q 

Let i\) = Re lS l h (R > 0) be a solution of the stationary Schrodinger equation 
(jZ|; assume for instance it corresponds to the iV-th energy level E]y. A direct 
calculation shows that the function R satisfies the equation 

— ~ S '( x )) R ( x ) = [ E ~ V(x)]R(x) (15) 

2m \ dx J 

and the corresponding Fermi function is thus 

g F (x,p) = (p- S'(x)f + 2m(E N - V(x)). (16) 

The area of the surface Q, bounded by the curve C : gp(%,p) = is given by 
the formula 

rx B 



1-xb 

Area(O) = 2 \/2m[E N - V(x)]dx. 

J X A 

In view of the quantization condition (|13p this means that we have 

Area(fi) = Nh + 2 / v / 2m[E - V(x)]dx. (17) 

J XA 
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Let us illustrate this on the harmonic oscillator with classical Hamilto- 
nian 

H = — (p 2 +m 2 uj 2 x 2 ). 
2m 

The ground energy is Eq = \ftw and the corresponding turning points are 
■iiy/h/muj hence 



ryjh/mu I — 

Area(fi) = Nh + 2V mHuj / Wl —x 2 dx; 

J-^fhf^L V "> 

the integral is easily evaluated, and one finds that 

Area(ft) = (N + \)h (18) 

and one thus recovers the fact that the action increases by jumps equal to 
h starting from the initial value \h. 

4.2 From Q to ip 

We now address the converse problem. Consider a smooth closed curve C 
(see Fig. 1) described by equations 

p = p + ( x ) , p = p~(x) (x A <x<x B ); (19) 

with p + (xa) = p~(xa), P + (xb) = p (xb) and p + (x) > 0, p + {x) < 0. 
Defining functions f(x) and g(x) through 



p + {x) = f(x) + yj-g{x) , p-(x) = f{x) - y/-g(x) (20) 
the curve C is given by the single equation 

(p-f(x)) 2 +g(x)=0 (21) 
and the area of the surface Q enclosed by C is 

Area(ft) = 2 f * v / -g(x)dx. (22) 

Jx A 

We are going to show that if C is quantized in the sense above, we can 
associate to it a solution ip(x) = Re lS ^l h of some stationary Schrodinger 
equation. Choose a number E > and define two functions S{x) and V(x) 
by 

f[x) = S'{x) , g(x) = 2m(V(x) - E) (23) 
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S(x) = j*f(x')dx' , V(x) = ^- + E (24) 



that is 



where xq is an arbitrary fixed number. Eqn. (|2ip for C is then 

^-(p-S , (x)) 2 + V(x)-E = 0. (25) 
2m 

Consider the differential operator 

2^(-<- S 'w) 2 + FW; 

it has a discrete spectrum consisting of positive eigenvalues. Let Eq be the 
smallest eigenvalue; the boundary of f2 must thus satisfy the quantization 
condition ([13]), that is 



/ y/2m[E - V(x)]dx = / y/2m\E Q ~ V{x)]dx + — 

JXA Jxa ^ 



IX A JXA 

This condition determines the value E. Let us now look for a function R(x) 
such that 

. , „ h 2 R"(x) . . 

V(x) -E = -j-f- 26 

that is 

- ^-R"(x) + (V(x) - E)R(x) = 0. (27) 
2m 

The function tp(x) = Re tS ^/ h is a solution of the equation 

2^ir ih i^- s ' {x) ) 2 +v{x) ^ (x)=0 

as is shown by a direct calculation using the identity (|27p ; ?/>(x) is, therefore, 
the wave- function we are looking for. Notice that the choice of the value xq 
in (|24|) is irrelevant because if we replace it with another value x' ^ xq it 
changes ip into tp' = e^^ip where 7 = f(x')dx' is a constant phase. 

5 Examples 

5.1 Squeezed states 

We consider (unnormalized) squeezed coherent states 

Mx) = e-^ a+ib > 2 (28) 
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where a and b are real and a > 0. The functions S(x) = — \bx 2 and 
R{x) = e - ax2 ' 2h satisfy 

s 'W = - te -^ = -5 + ^- (29) 

The Fermi function of ^a,b is thus the quadratic form 

gF(x,p) = (p + 6x) 2 + a 2 x 2 — ah. (30) 

Setting # = ( X ) we can rewrite formula (f3Q|) as 



g-p(x,p) = z T Mz — ah 
where M is the symmetric matrix 

Since det M = a 2 it follows that the surface O enclosed by the ellipse C : 
gF(x,p) = is i/t. 

A straightforward calculation shows that the matrix (|3ip factorizes as 



M = ^(° °) S ,32) 
where £ is the unimodular matrix 

^ = ( -l/2 ft -l/2) • ( 33 ) 

It turns out -and this is really a striking fact!- that Mp is closely related to 
the Wigner transform 

1 f°° 

Wijj a , h {z) = — J e-&y^ h {x + \y)^ h {x - \y)dy (34) 

of the state ip a ,b- I n f ac t (see e.g. de Gosson [5J, Littlejohn [15]). 

WlM*) = (Tr^-V^-i/^e-^fe/^ (35) 

where G is the matrix 

G = S T S=( a \ h ' /a b Jf). (36) 



ft/a V a y 



s 



It follows from Eqn. ([32]) that 

W^ a)b {z) = (7r^)- 1 / 2 (deta)- 1 / 2 e-°exp 
a 0' 



-\g ¥ {S-'D-^Sz) 

h 



(37) 



with D 

y a 

In particular, when n = 1 and ip a ,b{ x ) = e~ x ' 2 / 2h we have S~ 1 D~ 1 / 2 S = I 
and a = 1 and hence 

which was already noticed by Benenti and Strini [3] . 

5.2 Hermite functions 

The iV-th Hermite function is 

h N (x) = (^2 N N\)-^e-^l 2 H N {Q) , Q = x (38) 

where 



h 

H N (x) = (-l)^V |- e - (39) 

is the Af-th Hermite polynomial; the latter satisfies the second-order differ- 
ential equation 

H'x(x)-2xH' N (x)+2nH N (x)=0. (40) 

Here, S(x) = and a straightforward calculation using the relation (|40p 
yields 

+ (41) 

hence the Fermi function is here 

g F (x,p) = p 2 + m 2 uj 2 x 2 - (2N + l)muh. (42) 

The curve C : gp(x,p) = is again an ellipse, enclosing a surface Q with 
area 

Area(ft) = (N + \)h. (43) 

This example is not very instructive, because the Fermi operator gp is, 
up to the factor l/2m just H — {N + \)hjj where 



h 2 d 2 l 
1 — i 

2m dx 2 2 



H = — 7T H — moj 2 x 2 
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is the harmonic oscillator Hamiltonian, whose eigenstates are precisely the 
Hermite functions ()38|) . We leave it to the reader to verify that the same sit- 
uation occurs for every real function which is an eigenstate of some arbitrary 
operator 

h 2 d 2 



2m dx 2 
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